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We use the system-plus-reservoir approach to study the quantum dynamics of a bipartite con- 
tinuous variable system (two generic particles). We present an extension of the traditional model 
of a bath of oscillators which is capable of inducing an effective coupling between the two parts 
of the system depending on the choice made for the spectral density of the bath. The coupling is 
nonlinear in the system variables and an exponential dependence on these variables is imposed in 
order to guarantee the translational invariance of the model if the two particles are not subject to 
any external potential. The reduced density operator is obtained by the functional integral method. 
The dynamical susceptibility of the reservoir is modelled in order to introduce, besides a character- 
istic frequency, a characteristic length that determines if the effective interaction potential is strong 
enough to induce entanglement between the particles. Our model provides a criterion of distance 
for identifying in which cases a common environment can induce entanglement. Three regimes are 
found: the short distance regime, equivalent to a bilinear system-reservoir coupling, the long dis- 
tance regime in which the particles act like coupled to independent reservoirs and the intermediate 
regime suitable for the competition between decoherence and induced-entanglement. 
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INTRODUCTION 



The usual model of Brownian motion has been suc- 
cessfully used to describe general properties, classical or 
quantum mechanical, of dissipative systems with only 
one degree of freedom subject to arbitrary potentials 
[3, d, |j, @, 0|- Indeed, it has been extensively shown 
in the literature that, within the range of interest, other 
approaches to dealing with dissipative systems described 
by a single dynamical variable always furnish us with the 
same results as those obtained by the bath of oscillators 
with a properly chosen spectral function 
However, in a previous paper |l| was shown that the 
model of independent oscillators coupled bilinearly in co- 
ordinates to the system of interest is inappropriate to 
dealing with dissipative system in which two indepen- 
dent degrees of freedom are considered. There, a gener- 
alization of the usual model was developed extending the 
coupling to be nonlinear in the system variables and mod- 
ifying the spectral function to mimic the low frequency 
limit of the response function of an interacting oscilla- 
tors bath. The importance of understanding the dynam- 
ics of dissipative bipartite systems arises, for example, 
from quantum information and computation where it is 
mandatory to study the interplay between the decoher- 
ence and entanglement induced by the environment. 

It is our intention in this paper to investigate the quan- 
tum behavior of a system composed by two independent 
particles immersed in a common environment. We use 
the path integral approach to find the reduced density 
operator, in coordinate representation, of the two parti- 
cle system. The propagator that controls the temporal 
evolution of the system is written, as usual, in terms of 
the influence functional that contains the effects of the 



bath. In order to guarantee that each particle behave 
as a Brownian particle, if the other is absent, the in- 
fluence functional can be expressed as a triple product 
where two factors are equivalent to the result one would 
have obtained by coupling each particle separately to a 
bath of noninteracting harmonic oscillators in the usual 
way and the third factor introduces and effective inter- 
action between the particles mediated by the reservoir. 
The latter behaves as a source of quantum correlations 
and allows us to find, for example, entangled states for 
systems initially prepared in a separable state. 

As an application we calculate the temporal evolu- 
tion of the entanglement (logarithmic negativity) for a 
Gaussian two mode state. The results show a competi- 
tion between decoherence and induced entanglement that 
mainly depends on the temperature and the average dis- 
tance between the particles. 



MODEL AND EXACT DENSITY OPERATOR 

In a previous paper we presented a generalized 
system-reservoir model suitable for the study of two 
uncoupled particles interacting with a common bath. 
There, it was demonstrated that the traditional bilin- 
ear coupling is not appropriate to handle the effective 
interaction mediated by the bath and the solution was to 
introduce a nonlinear coupling in the system variables. 
The principal features of the model will be reviewed in 
this section. 



Model 



The Hamiltonian for the complete system is given by 
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One Particle Density Operator 



Hs + Hr + Hi. 



(1) 



Hs is the Hamiltonian of the system of interest, Hr is the 
Hamiltonian of the reservoir, which will be described as a 
symmetrized collection of independent harmonic modes 
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and Hi is the interaction Hamiltonian which can be writ- 
ten in two equivalent forms 0, 0] ; a coordinate- velocity 
coupling or a coordinate-coordinate coupling plus a new 
quadratic term which is necessary to preserve the transla- 
tional invariance of |T|) when the system of interest is not 
acted by any external force. The two forms are linked by 
a simple canonical transformation (P k — > m k u> k R k and 

In our generalization of the previous model, the sys- 
tem of interest with a single degree of freedom will be 
represented by the free particle Hamiltonian 
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For the coupling term we assume the interaction 
Hamiltonian 

1 N 

Hi=-J2 (C~k(x)R k + C k {x)R^ k ) 
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In order to represent the effect of a local interaction of 
the particle with a spatially homogeneous environment 
we choose 



C k (x) 



ikx 



(5) 



With this choice it is straightforward to show that the en- 
tire system is translationally invariant when there is not 
an external potential [l| . As a consequence of ([5]) the po- 
tential renormalization in Q is a constant and therefore 
does not contribute to the particle dynamics. Actually, 
this reflects the translational invariance of whole system. 
A coupling like (0 appears, for example, when one deals 
with the interaction of a particle with the density opera- 
tor of a fermionic bath 
problem 

The classical equations of motion and the implications 
of the coupling §5§ were studied in 1] for one and two 
particles. There the nonlocal influence of the bath ap- 
pears only when a second particle is present in the en- 
vironment. Here we show that the nonlocal effects are 
also important for a single particle when the quantum 
behavior is studied but the local form [2] is recovered in 
a suitable limit. 
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in the traditional polaron 



In this section we use the functional integral method 
developed by Feynman and Vernon [5( in order to find 
the reduced density operator of the system of interest. 
The time evolution of the total density operator is given 
by 

p(t) = exp (-iHt/h) p(0) exp (iHt/Ti) , (6) 
and in coordinate representation we have 



(x,TL\p(t)\y,Q) 



dx'dy'dR'dQ'K(x, R, t; x , R , 0) 



*(x',R'\p(0)\y',Q')K*(y,Q,t;y',Q',0), 

(7) 



where R is a iV-dimensional vector representing each de- 
gree of freedom of the bath and K is the coordinate rep- 
resentation of the time evolution operator. The reduced 
density operator emerges when we eliminate the bath 
variables. In the coordinate representation this proce- 
dure corresponds to doing R =Q in (|7|) and to integrate 
over all possible values of R. The reduced density op- 
erator at time t depends on the total density operator 
at t = and for simplicity we assume a separable initial 
condition 



p(0) = p(0) Pfl (0), 



(8) 



where p(0) is the initial particle density operator and 
Ar(0) = Z^ 1 exp(—(3Hn) is the bath density operator, 
which is assumed to be in thermal equilibrium, before 
the perturbation is switched on. With this assumptions 
we can write 



p{x,y,t) 



dx'dy'J(x, y, t; x' , y' , 0)p(x', y', 0), (9) 



where J(x, y, t; x', y', 0) is the superpropagator, which 
controls the time evolution of the reduced density op- 
erator and can be written as a path integral 

x y 

J(x,y,t;x',y',0)= [ Vx{t') [ Vy(t')* 



exp - {So W)] - So [y(t')}} T [x(t% y(t')] . (10) 

So is the action of the isolated particle and J- [x(t'), y(t')] 
is a functional of the particle trajectory, which contains 
all the bath information [5| . The influence functional has 
a simple formal expression 

F[x(t'),y(t')] = Tyr ( P RU RI [y(t')pRilx(t')}) , (11) 

where Um[x(t')] is the unitary time evolution operator 
of the reservoir subjected to the influence of the system 
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which evolves through the Hamiltonian H R j = H R + 
Hj[x(t')]. This means that a given trajectory x(t') of 
the system for < t' < t acts as a forcing term to the 
environment. The time evolution of the operator U R i[x] 
is described by the Schrodinger equation 

iKjU R i(t) = H RI (t)U RI (t), (12) 

with the initial condition U R j(0) = 1 and has the formal 
solution 

U RI {t) =Te- i ti dt ' HRl( - t ">/ n , (13) 

where T is the time ordering operator. Switching to the 
interaction picture [l3| this result can be written as 

U R i{t) = e -^ R t/h Te - l f ( t > dt'H I [ x (t')]/n^ (w) 



where Hi[x(t')] = e iHRt / n Hi[x(t')]e- iHRt / n . In- 
serting CH} mto CP we have T [x(t'), y(t')] = 
Tr R {pRTetft^ZthW^Te-ifi^^VwA . 

Since the reservoir have many degrees of freedom we 
can assume that the particle induces a weak perturbation 
in the environment and expand the chronological product 
to second order in Hj [Icj]. 

Te ^^dt'H l[:c (t')]/n ra i _ * f dt'HMt')] 

Jo 

t t' 

-~2 I d£ I dsHi[x(t')]Hi[x(s)]. (15) 
fi Jo Jo 

With these assumptions and tracing the reservoir vari- 
ables we find 



F[x(t'),y(t')} » l-^j\t'J^ ds[{M^nM<s)}) + {Mv{s)]Hv{n) 
(HjW^HMs)]) - (Hi[y( S )]Hi[x(t')})) . 
Once we consider only terms to second order in Hi[x(t')] the influence functional above is equivalent to 

T[x(f),y{t')] = eW i~J*dt' jf* ds((ffj[x(0]#/[»(»)]) +(#/[»(»)]#/ ly^)]) 

- (HMt'MMs)}) - (Htiy^HMt')}))} ■ (16) 

I 



Now we calculate the averages in (fT6|) considering the 
coupling |5|) and invoking the reservoir translational in- 
variancc. Since we have (Rk(t')Rk'(s)) = unless k' — 
—k, one is allowed to write 

(HMt')}Hi[x( S )}) = i£{C_ fc [s(0]C fc [z(«)] 

+C , fe [ a3 (t / )]C-fc[^(s)]} (R k (t')R- k (s)) . (17) 
The expression (Rk{t')R-k{s)) can be written in terms 



of the bath linear response through the fluctuation- 
dissipation theorem. We call for simplicity ctk(t' — s) = 
(Rk(t')R- k {s)) and 

ft, f°° e -iu(t'-s) 

a k {t' -s) = - J ^ du>Imxk(u) 1 _ e _ uH/3 , (18) 

where Xfe(w) is the dynamical susceptibility of the reser- 
voir. So, the averages in (fl6|) have the final form 



I 

(fli[x(f)]flj[a:(fl)]) = i^{C_ fc [x(0]Cfc[x(s)] + C fc [a!(0]C'-fc[ir(»)]}a fc (* / -«) ) (19) 

and similar ones for the other terms. Putting all those together we find the following nonlinear influence functional 

!F[x(t'),y(t')] = expj — / di' / [w^K-fc (cos k{x(t') — x(s)] — cos k[y(t') — ctfcft' — s) 

I ft Jo Jo k 

+ n k ^ k (cos k[y{t') - y{s)} - cos k[y{s) - x{t')])a\{t' - a)] \. (20) 
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The function ctk(t' — s) is related to the bath linear 
response by the Fourier transform (|18p . At this point it 
is necessary to model the bath dynamical susceptibility 
Xfe(w), since it is not our intention to completely describe 
the microscopic details of the bath. We can assume the 
imaginary part of the bath dynamical susceptibility has 
the form 



Imxfe(w) = f{k)uj6{uj - O). 



(21) 



Here we introduce a high frequency cutoff f2 as the char- 
acteristic frequency of the bath. The Markov dynamics 
is achieved when we take the limit — > oo, and the func- 
tion f(k) responds for the nonlocal influence of the bath. 
This approximation is equivalent to replace the free oscil- 
lator response function by the low frequency limit of the 
damped oscillator response function [l| and allow us to 
separate the characteristic time and length scales of the 
reservoir. A functional dependence like (f2"Tj) for the dy- 
namical response of the bath has been employed in Refs. 
ToL 11 1 for fermionic environments. 
The parity of Imxfc(w) allow us to write the correlation 



function otk(t' — s) as 

a k (t'-s) = a[ R) (t' - s) + ia™ (f - s) , (22) 
with the real and imaginary parts defined respectively as 



( R ) 14.1 \ 

n 



7T 



dwlmxfc(w) cosuj{t' — s) coth(7i/3u;/2) (23) 



and 



du)Im.Xk{u) sinu;(i' — s).(24) 



The functional (|20|) can also be written in terms of these 
real and imaginary parts. With prescription (|2ip we can 
evaluate the frequency integrals in the imaginary part 
(f2~4")l taking the limit f2 — > oo and considering that the 
coupling was switched on at t = + . Following this pro- 
cedure we have a functional with a Markovian imaginary 
part and a explicit nonlinear dependence, 



F[x{t'), y(t')] = exp j JL ^ KfcK „,/(fc)fc jf dt 1 sin k[y{t') - x{t')] {x{t') + y{t')) 

— "2 / dt' f ds KkK-k [cos k[x{t') — x{s)\ — cos k[y(t') — x(s)] 
h Jo Jo k 

+ cos k[y(t') - y(s)} - cos k[y(s) - x(t')}} a[ R) (f - s) 



r 



This functional can be approximated assuming that the example, in fermionic environments this length is related 
most important trajectories (x(t'),y(t')) are confined to the Fermi wave number fcj? U|- In this approxi- 



within a region small compared to a characteristic length mation, we have k{y{t') — x(t')) <C 1 and the functional 
k^ 1 introduced in (j2"Tj) through the function f(k). For can be written as 



T[x{t'\ y(t')} = exp |l| J* dt'(y(t') - x(t')) (±(f ) + y(t')) 
v f l , / /"*' , f°° , , Chfr 

dt ds duiui coth — 

tin Jo Jo Jo V 2 



I" dt' [ ds [ dujuj coth fO^. \ ( x (t') - y (t')) cos oj(t' - s)(x(s) - y(s)) } , (25) 



r 



where we have identified 77 = J^k k 2 KkK-kf{k). Notice the result obtained by coupling the particle of interest 

that, with this modification, we obtain a relation between bilinearly to a bath of nonintcracting harmonic oscillators 

the damping constant and some microscopic parameters with the spectral function J(u>) =770; 
of the oscillator bath. The functional J25[) coincides with 
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Two Particle Density Operator 

Now we are going to study the dynamics of a system 
with two degrees of freedom immersed in a dissipative 
environment. In this case the Lagrangian of the system 
of interest is 



1 



~-Mx\ 



1 



Mil 



2 1 2 
and the coupling term 

Li = -\Y. [(CU(zi) + C- k (x 2 )) R k 

k 

+ (C k (xi) + C k {x 2 )) R_ k ] . 



(26) 



(27) 



Notice that we have not included any counter-term in 
(|27[) since our system is manifestly translationally invari- 
ant. The time evolution of the reduced density operator 
for the two particle system is given by 

/5(x,y,t) = JJd 2 xd 2 yJ{x,y,t : x',y',0)p(x',y',0) 7 

(28) 



and the superpropagator for this case is 
x y 
J(x,y,t; X ' ) y' ) 0) = J Px(t') J Dy(i') 



exp 1 {So [x(f )] -So W)]} F W)Mt')} , (29) 

where we have defined the vectors x(£') = (x\(t') , x 2 {t')) 
and y(t') = (yi(t'),y 2 {t')). So is the action of the iso- 
lated two particle system and F is the Feynman- Vernon 
influence functional, which in the operator form can be 
written as 



F 



[x(i'),y(0] = Tr* (pRUhbWpmW)]) , (30) 



where Uni[x(t')] is the unitary time evolution operator 
as in (flUj) . The procedure to calculate the influence func- 
tional for the two particle case is completely equivalent 



to the one particle case and we can directly generalize the 
result, keeping only terms to second order in Hj[x(t')], 
as 



ynx(/'l.yU')]=rxi><!~ 



t r-t 

dt' I ds 
Jo 



x ^JJ z [x(f )]ffj[x(*)]) + {HfrWHjfyV)]) (31) 
>,[y(f)]tf/[x(*)]\ -/^ I [y( fl )]fl>[x(t')]\)V (32) 



Again, we take the average over the equilibrium density 
matrix of the bath and using the fluctuation dissipation 
theorem and the correlation function defined in (J22J), we 
can write the two particle influence functional as 



F[x(f),y(f)] = 



expl j / dt ' / dsK[x, y, t', s]al ' (t 1 - s) 

\ fi Jo Jo 



t rt 

dt' I dsK[x,y,t',s]a k I} (t' - s) 



(33) 



fl JO JO y 

where we have defined the nonlinear kernel 
K[x,y,t', s] = 

2 

^KfcK_ fe ^2 [cosklxitf) -Xj(s)] + cos k[yi(t') - yj(s) 

k i,j=l 



cosk[yi(t') — Xj(s)] — cos k[xi(t') — yj(i 



(34) 



The imaginary part of (|33|) can be reduced using the 
response function (|2ip in a form completely analogous 
to the one particle case. The resulting imaginary part of 
the exponent of the influence functional can be written, 
in a compact form, as 



1 



21) 



29. 



— dt'L[x,y,t'} + — dt'V[x,y,t% (35) 



Hit 



'0 jo 
where we have defined the instantaneous kernels 



L[x,y,t'] =J2 K kK- k kf{k)( ^ WRk[xi(f)- yj (t')] [i<(0 +%(*')] 
k Vi=l 

+ smkixtit') - x 2 {t')] [x 2 {t r ) - ii(i')] - sin k[ yi (t') - y 2 (t')] [y 2 (t') - yi {t')] 



(36) 



and 



V[x,y,t'} = «*«-*/(*) (cosfeMO ~ *a(0] ~ cos k[ yi (t') - y 2 (t')]) . 



(37) 
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The well-known effects of dissipation and diffusion are 
included in the kernel L[x,y,i']. Note that this term 
contains direct and indirect influence of dissipation . On 
the other hand, the function V[x, y, t'\ introduces a com- 
pletely new effect that can be interpreted as an effective 
interaction mediated by the environment. The functional 
at this point can be written as 

F[x(f),y(f)] = 

exp < 2 / dt' I dsK[x.,y,t' , s]a^ (t' — s) 

[ ft Jo Jo 

-1- f t dt'L[x,y,t'} + i ^ f dt'V[x,y,t>}\ . (38) 



2h 



o 



Hit 



The functional above is still too complex and addi- 
tional assumptions are necessary. Since we are con- 
cerned only with the effective terms mediated by the 
reservoir, we can assume that each particle trajectory is 
localized within a region restricted by the characteristic 
length of the reservoir. That leads to the approximation 
k[xi(t') — Ui(s)] <C 1 and taking only the terms up to sec- 
ond order in k[xi(t') — yi(s)] the functional can be written 
as 

F [x(f ), y(t')} = T [si(f ), W (f )] T [X2(t% y 2 (t')} 

exp i — \ [ d£ [ dsK[yi, y, t' , s]a[ R) (t' - s) 
{ ft Jo Jo 

-± f t d t'L[ X ,y,t'} + ^ [* dt'V[x,y,t'}\ , (39) 



2ft 



ftlT 



where K and L are defined as in (I34|) and (|36|) . respec- 
tively, but excluding the terms with i = j. The fac- 
tor T[xi{t'),yi{t')\ was defined in (125)) and corresponds 
to the influence functional within the local approxima- 
tion of the Brownian particle dynamics. The functional 
above is formed by three clearly distinct parts, the first 
two correspond to the direct influence of the bath over 
each particle separately and the third part contains the 
interaction between the particles mediated by the bath. 



Towards a better understanding of the interaction fac- 
tor in (|3"9")) it is important to study the short and long 
distance limits, which are obviously defined in relation 
to the characteristic length of the reservoir k^ 1 . For the 
first we assume k[xi{t') — Xj(s)] <C 1 for all particle tra- 
jectories and for all times, which means that the particles 
are very close and the dynamics is essentially local. In 
this approximation the superpropagator reduces to 



x y , 

J = e X J^ft(X,Y)) J j T>Ktf)Ttytf)expUSo[xtf)]-S \ytf)]-2 J fcitfi - Vi±i)dt' - | J (x 2 y 2 - y 2 x 2 )dt' 



- 77 / dt ' [ x ih - 2/2*i + x 2 yi - y%x 2 ] -—[ dt' (xi - x 2 f - (yi - y 2 ) 2 

2 Jo T Jo L 

xexpj-^jf dt'jf dsJ2 [x i (t , )-y i (t , )}K(t'-s)[x j (s)-y j (s)} 



r 



where we have introduced the temperature dependent 
kernel 

K(t'-s)= / du uj coth(ft/3uj / 2) cos uj(t' - s) (40) 



and the function 



h,(X,Y) = ^[( Xl +x 2 f-(x' 1+ x' 2 ) 2 

2 i /„./ i „/\21 



-{yi + V2) 2 + {y' 1 +y' 2 f 



(41) 



Now, for simplicity, we can use two successive changes of 
variables. The first is defined as qt(t) — ( Xi (t) + y i {t))/2 1 
£i(t) — Xi{t) — yi(t) and can be interpreted as the cen- 
ter and width of the wave packets respectively. The sec- 
ond introduces the center of mass and relative coordinate 
variables defined respectively as r(t) — (qi(t) + q 2 (t))/2, 
u(t) = qi(t) — q 2 (t) and the auxiliar variables x(t) — 
(£i(i)+£ 2 (t))/2, v(t) =£i (*)-&(*)• With these replace- 
ments the local version of the superpropagator reads 
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J = exp \ —h(r,x) 



r X u v 



■Dr(t')Vx(t')T>u(t')T>v(t') exp j f dt' / <ls K f /' - ,j \ ( /' i \ ( , ) 



r' x' u ' v ' 
rt 



JO 



exp-<V[r, x ,u,v]-2T} / dt'[ X (t' )r(t') - r(t ')*(*')] 



dt'u(t')v(t')dt' 



(42) 



with 

Z[r,x,u,v]= [ M(2rx + uv/2)dt' (43) 
Jo 

and 

fc(r,x) =2J7(rx-r'x'). (44) 

From the local functional (|42|) we can point out the fol- 
lowing features. Firstly, the stationary trajectory of r(t), 
deduced from the action functional in the propagator's 
exponent, describes the dissipative dynamics of a par- 
ticle with mass 2M, that is the center of mass of the 
two Brownian particle system. The second, and more 
remarkable result is that, the relative coordinate u(t) in 



this approximation describes a dissipationless dynamics 
of a variable acted upon by an effective force which in- 
duces an interaction between the individual parts of the 
system of interest. This effect had previously been no- 
ticed in a classical approach and coincides with that of 
other works in which the quantum evolution was stud- 
ied using the master equation approach and a bilinear 
system-bath coupling 15J, |l6[ . 

In the long distance limit we suppose that each par- 
ticle is restricted to move within separated regions. In 
this case the distance L between the regions is consid- 
ered bigger than the characteristic length of the reservoir 
and it is possible to approximate k[xi(t') — yi(s)] ^> 1. 
When this condition is satisfied the interaction terms are 
negligible and the superpropagator reads 



T U X V 



J = 



Vr{t')Vu{t')Vx{t')Vv{t') exp { - / dt' M 



r' u' x' v ' 



exp 



iV f* dt , fu(t')v(t') 



+ 2r(t')x(t') - 



Tin 



u(t')v(t') 



2r(t')x(t') 



'dt' f ds f^M + 2 U .s A (/')) Kit' -si 
o Jo 



(45) 



It is obvious that in the long distance limit the dynamics 
is completely equivalent to the case of two uncoupled 
particles interacting with two independent but identical 
environments. 

Now, in order to study an approach suitable for in- 
terparticle distances between the two above-mentioned 
limits, it is necessary to introduce new assumptions. We 
consider the variables insensitive to the interparticle dis- 
tance negligible in relation to the characteristic length 
/cq -1 , i. e. the variables v and x j which are related to the 



width of the wave packets, can be kept as very small when 
compared with k and the approximations cos [fcx(t')] ~ 1, 
sm[k X (t')] w kx(t'), sin[kv(t')/2] w kv(t')/2 appear to 
be adequate ones. On the other hand, the nonlinear 
terms involving the variables r and u will be replaced by a 
phenomenological function parametrized by the relation 
between the characteristic length of the bath and the av- 
erage distance between the physically allowed regions for 
the particle trajectories. In this way the superpropagator 
can be written as 
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J = exp<j -h(r,x,u,v) 



r U X V 



Vr(t')Vu(t')V x (t')Vv(t') 



r' u' x' v ' 



xexp^A// dt' (2rx + uv/2) + 1 l f dt' (4x(t>(t') + v(t')u(t') - u(t')v(t') - 4r(f')x(*')) 



dt' (D(k L)v(t')u(t') + 4D(k L)r(t') X (t')) 



2 V n f* 



dt' D(k L)u(t')v(t') 



x exp 



J* dt' jT dsK(t' - s) 1 2 (1 + D(A„i)) X(*')x(s) + 5 (1 - £>(*Qi)) «(*>(«)} 



(46) 



r 



where we redefine 



h(r, x, u, v) = — (4r x — 4r% + it v — uv) 



n 



D(k L) (u'v' - uv) , 



and introduce the function 

D(k L) = cxp (-k L) , 



(47) 



(48) 



with the parameter L being a definition of the average 
length between the regions where the particles move. Ac- 
tually one could use it as the initial distance between the 
centers of the packets once one can make sure they will 
not overlap in the long run. 

The functional integrals in (|4*6"|) can be evaluated di- 
rectly since all the variables appear to second order only. 
The result of this integration is 



J = Af { t)e Xp Uh(r,x,u,v)~-£^ M 



M 



7+ (r - r) (xe- 1+t/2 - x'e 7+t/2 ) - 7 (rx ~ H) [*y + /2] 



4 sinh 



(7_ — 7) (uv — uv) sinh 
1 



y^ 2 - 4lj? (u've-^ 2 + w'e^"*/ 2 ) 

— 4o> 2 (uv + u'v') cosh ^ \Jl- — 4wf I 



x exp--|(l+£)(fcoi)) (A x (t) x 2 + B x (t)xx' + C x (t) X ' 2 ) + (1 - (A,ft)w 2 + B„ft)W + C„ftK 2 ) j. (49) 

In this expression we use the definitions Uq = -fp 2 , 7± = 7(1 ± D[fcoi]), <jj\ = u>QD(koL) and the time dependent 



functions are 



A x(t) = - / dr / dsK(T-s) 



Jo 



,-7+* 



sinh [S7+/2] sinh [t7+/2] 
sinh[t7 + /2] 2 



,7+(s+r)/2 



(50) 



B x( t ) = — I dT I dsK(r-s)e 



Jo 



,-7+*/2 



sinh [r 7+ /2] sinh [ft - s)7+/2] 
sinh [i7+/2] 2 



l7+ (s+r)/2 



(51) 



C x (t) = - / dr / dsK(r-s 



Jo 



sinh [ft - s)7 + /2] sinh [(t - r) 7+ /2] 7+(s+r)/2 



sinh [*7+/2] 



(52) 



AM) = — / dr / dsK( T -s)e 
4tt Jo 7o 



^ /-j.w. -N--y-t L 2 V LL_ L 2 V- M 7 _( a +r)/2 



sinh I \/7- — 4w 2 sinh § a/7 2 . — 4w 2 



sinh 



(53) 
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B v (t) 



2tt 



dr 



t sinh 
dsK{r - s)e-^- t/2 



sinh 



(ts) 
2 



4u; t 2 



sinh 



7 2 -4 W 2 



o7 -(s+r)/2 



(54) 



= — / dr / <fcKTr-«)- 
4tt Jo Jo 



sinh 



4u, t 2 



sinh 



4 W 2 



sinh 
I 



7 2 - 4w t 2 



o7 -(s+r)/2 



(55) 



jV(t) is a time dependent coefficient resulting from the 
fluctuation around the classical path and will be deter- 
mined imposing the normalization condition to the final 
reduced density operator. With the propagator (|49[) we 
have all the tools needed for studying the time evolution 
of the reduced density matrix which we do in the next 
section with a particular example. 



EVOLUTION OF A TWO PARTICLE STATE 

Now, we want to study the consequences of the ef- 
fective interaction induced by the reservoir in the time 
evolution of a system formed by two particles, without 
any direct interaction between them. Our atention will 
be focused on the entanglement dynamics of a bipartite 
system coupled to an oscillator bath and, for clarity, a 
brief introduction about how to measure entanglement 
in Gaussian states is necessary. 

A Gaussian state can be described using the charac- 
teristic function, which in general can be written as 



H*(X) = cx V \ -ixAX 5 



(56) 



where X = (i>\, Ai, i>2, A2) and A is the covariance matrix 



with elements A 



(XiXj 



XjXi 



/2. The character- 



istic function is related to the density matrix through the 
transform 



exp (— iXx./Ti) (x 



P 



d x. 



(57) 

where A = (Ai, A2) and u — (v\,V2) are a sort of momen- 
tum and position variables respectively. The covariance 
matrix contains all the useful information of the Gaus- 
sian state and determines it completely. The covariance 
matrix shall be conveniently written in terms of 2 x 2 
block submatrices as 



A = 



A C 
C T B 



(58) 



A very common quantifier of entaglement is the logarith- 
mic negativity E^f defined by 



Ejsf(p) = max [0, — ln2cr_] 



(59) 



where a± are the symplectic eigenvalues of the partial 
transposed density matrix that can be written using the 
local symplectic invariants detA, det/1, det-B, detC as 



1 

V2 



D A ± 



4 det A 



(60) 



where D\ = det A + det B — 2 det C. This expression 
quantifies directly the violation of the PPT (positive par- 
tial transposition) , a necessary and sufficient condition of 
separability. 

We choose as initial state a Gaussian density matrix of 
two modes caracterized by the squeeze parameter z. The 
initial density matrix written in the variables r, u, v, Xi 
is 



p(r', X >'y,0) 



2tto- 



■ exp 



-2z 



4ct 2 



(4/ 2 + X ' 2 ) 



„2.- 



4ct 2 



{u' 2 +v' 2 /A) 



(61) 



The density matrix at t > is obtained from the expres- 
sion PB"]) , where the propagator J and the initial density 
matrix are given respectively by (|49p and (|6ip . with the 
suitable change of variables. The integrals indicated in 
(|28| can be solved exactly for all the limits mentioned 
above since the involved variables appear only up to sec- 
ond order. After the evolution of those integrals the co- 
variance matrix is easily computed by the Fourier trans- 
form ([57)) and the results of l|59p are presented for various 
parameter combinations in the figures below. 

First, lets us observe the evolution of the logarithmic 
negativity when the two parts of the system of interest 
are in the local (short distance) regime. In figure fl} we 
have the logarithmic negativity as a function of time for 
low temperatures and three values of the squeeze param- 
eter. The dependence on the initial entanglement is evi- 
dent in this figure. For a initial separable system (z = 0) 
the environment rapidly induces entanglement and on the 
average maintains a value higher than that for systems 
initially entangled in which the decoherence process dras- 
tically reduces the logarithmic negativity. The oscillatory 
behavior is due to the effective interaction potential in- 
duced by the environment [H and it is remarkable that in 
all the cases plotted in figure |T]) there exist a remanent 
finite entanglement. 
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Figure 1: logarithmic negativity for close particles as a func- 
tion of time, the parameters are specified at the top of the 
figure. 



Figure 3: logarithmic negativity for close particles and finite 
temperature as a function of time. The parameters are spec- 
ified at the top of the figure . 



y=1.0; kT=10.; Jt L=0.; H=10 



z=0 

z=0.5 

z=1.0 
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Figure 2: logarithmic negativity for close particles as a func- 
tion of time. The parameters are specified at the top of the 
figure . 



Figure @ show us that the remanent entanglement is 
an effect associated with low temperatures and when the 
temperature is raised the initial entanglement always go 
to zero. Again in this case, the initially separable state is 
the most robust and maintains the induced entanglement 
for longer times. 

When the particles are very close together the most 
important influence of the reservoir is the effective inter- 
action induced between them and the decoherence effect 
only becomes important when the temperature rises. In 
figure © the logarithmic negativity is plotted for three 
values of the dissipative constant 7. The average amount 
of entanglement is higher for higher values of 7, i.e. 
stronger interaction with the reservoir. Results similar 
to this, but following other formalism and always con- 
sidering a bilinear system-reservoir coupling, have been 
presented by other autors 15 1. 



Up to this point the results showed correspond to bi- 
partite system for which it is possible to define an average 
distance L between the two parts and this distance is very 
small when compared with the characteristic length of 
the reservoir, fc^ 1 . This approximation assumes that the 
relative coordinate is completely free from any dissipa- 
tive influence. In a more realistic situation a dissipative 
dynamics for the center of mass and relative motions is 
expected and it is plausible to find a monotonically de- 
creasing remanent entanglement when the distance be- 
tween the particles is increased. In figure (U]) we can 
observe the influence of the distance between the parti- 
cles on the evolution of entanglement. It seems there is a 
competition between the process of entanglement induc- 
tion and decoherence. In particular, when the distance 
increases the average entanglement is reduced and the 
oscillations become damped. In the curve for koL = 2.0 
a delay in establishing entanglement can be seen. In- 
creasing the distance L reduces the induced interaction 
and increases the delay time. In the limit k^L — > 00 
the delay time goes to infinity and for finite distances an 
asymptotically entangled state when the temperature is 
low enough is observed. 



If the inital state is entangled the delay time is a revival 
time. The initial entanglement is lost by decoherence due 
to the action of the reservoir and after that the system is 
again entangled by the induced interaction. This effect is 
only appreciable at low temperatures. Then, it is possible 
to tunc the induced entanglement process varying the 
parameters and, in particular, controlling the distance 
between the subsystems as shown in figure (|S|). 
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r =1.0; z=0.; kT=10" 4 ; H=10 

E N 

2.0 r 

k L=0. 



yt L=0.5 




Figure 5: Influence of the distance on the revival of quantum 
correlations 

CONCLUSIONS 

In this paper we presented a system-plus-reservoir 
model with a nonlinear coupling in the system coordi- 
nates which, in the adequate limit, allows us to repro- 
duce the phenomenological results known for the quan- 
tum dynamics of a Brownian particle. The model was 
extended for dealing with a system of two uncoupled par- 
ticles immersed in a common environment. The choice 
made for the behavior of the bath response function, in 
order to reproduce the dynamics of the quantum brow- 
nian motion for each particle when isolated, naturally 
generates an effective interaction which depends on the 
average distance between the particles and the dissipa- 
tion constant 7, as can be seen from the exact density 
operator. When the particles are very close together the 
interaction is basically harmonic and the amplitude of the 



induced quantum correlations oscillates without decaying 
when the temperature is low. In this approximation the 
dynamics is local and the system preserves the induced 
entanglement. It is worth noticing that when the initial 
state is entangled, the system evolution is less sensitive 
to the induced correlations and more sensitive to the de- 
coherence effects. For interparticle distances comparable 
with the characteristic length k^ 1 of the reservoir other 
phenomena appear, in particular we have the possibil- 
ity controlling the time interval between the total lack of 
initial entanglement and the induced entanglement. 

The results presented here may be relevant for a better 
understanding of the bipartite systems with regard to the 
theory of quantum information and the extension to deal 
with many Brownian particles might contribute to the 
study of quantum correlations of multipartite systems, 
and of the behavior of complex many particles systems 
in general. 
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